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A connection between quantum antiferromagnetism and high Tc superconductivity is theoreti- 
cally investigated by analyzing the t-J model and its relationships to the Gutzwiller-projected BCS 
Hamiltonian. After numerical corroboration via exact diagonalization, it is analytically shown that 
the ground state of the t-J model at half filling (i.e., the 2D antiferromagnetic Heisenberg model) 
is entirely equivalent to the ground state of the Gutzwiller-projected BCS Hamiltonian with strong 
pairing. Combined with the high wavefunction overlap between the ground states of the t-J model 
and the projected BCS Hamiltonian at moderate doping, this equivalence provides strong support 
for the existence of superconductivity in the t-J model. The relationship between the ground state 
of the projected BCS Hamiltonian and Anderson's resonating valence bond state (i.e., the projected 
BCS ground state) is discussed. 

PACS numbers: 



I. INTRODUCTION 

There are many reasons as to why high Tc supercon- 
ductivity has attracted so much attention. Aside from 
the obvious (but very important) prospects for techno- 
logical applications of high Tc materials, a very salient 
reason is the mysterious nature of the pairing mechanism. 
Since the bare electron-electron interaction is strongly re- 
pulsive, it seems that there are two ingredients essential 
to electron pairing: (i) how to overcome the Coulomb 
repulsion and (ii) how to generate an attraction between 
quasiparticles. 

In "low" Tc superconductivity described well by the 
standard BCS theory, the pairing mechanism is roughly 
as follows: Electrons form a Fermi liquid in which 
the strong Coulomb interaction is renormalized into 
a screened interaction between dressed quasiparticles 
which interact very weakly; with the strong repulsion 
gone, quasiparticles can form pairs through the exchange 
of phonons, which gives rise to a time-delayed attraction 
between quasiparticles. This intuition, while valid for 
standard BCS superconductors, is not correct for high 
Tc materials. 

First, in contrast to low Tc superconductors which are 
metallic in the normal state, cuprates arc insulators at 
low doping, and thus it is not a priori clear whether high 
Tc superconductivity in cuprates has anything to do with 
the Landau- Fermi liquid description of the normal state 
in metallic systems. In fact, this suspicion is reinforced 
by many non-Fermi liquid behaviors of cuprates includ- 
ing pseudogap phenomena 1 and stripe excitations^. In 
other words, it is not clear how the strong Coulomb repul- 
sion can be overcome by standard quasiparticle screening 
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mechanisms. 

Second, there is a wealth of evidence suggesting that 
pairing might be induced by a source other than the 
phonon exchange. While certainly the high energy scale 
of Tc in the cuprates itself is difficult to explain in terms 
of the phonon exchange mechanism, probably one of the 
most persuasive pieces of evidence for a non-phonon pair- 
ing mechanism is the destruction of superconductivity (or 
the strong suppression of Tc) when even a small concen- 
tration of Cu atoms (in the Copper oxide plane) are re- 
placed by non-magnetic impurities such as Zn 3,4,5 . This 
is very suggestive of pairing with magnetic origin. Also, 
very important in this context is the d-wave symmetry 
of the gap functionSiS*^, which is explained most nat- 
urally in terms of pairing due to magnetic interactions. 
Therefore, it seems reasonable to assume that magnetic 
interactions have something to do with high Tc super- 
conductivity. However, the precise relationship between 
magnetism and high Tc superconductivity is not very 
well understood; it is not clear exactly how the magnetic 
interaction generates an attraction between electrons, es- 
pecially at low doping. Also, from a fundamental point 
of view, this attraction is very puzzling since, after the 
phonon-exchange mechanism (i.e., interaction between 
electrons and ions) is ruled out, it is hard to appreci- 
ate how the remaining microscopic interaction could be 
substantially different from the completely repulsive bare 
electron-electron interaction. 

Despite the seemingly unrelated nature of the above 
two ingredients, they are, in fact, very closely connected. 
Specifically, it can be shown that the magnetic (antifer- 
romagnetic, to be precise) interaction is a natural conse- 
quence of the strong repulsive interaction between elec- 
trons at low doping, a process known as super-exchange. 
The precise mathematical derivation of super-exchange 
can be carried out in the framework of model Hamilto- 
nians. To this end, consider the Hubbard model with 
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repulsive on-site interaction U: 

H H ub = - tS ^( c \a c j<? +H.c.) +U ^2n i] n ii (1) 

where a =],[ is the spin index, and indicates that i 
and j are nearest neighbors. All models studied through- 
out this paper are defined on the square lattice. It has 
been showni2iii*i2ii^that, exactly at half filling, the Hub- 
bard model becomes identical to the antiferromagnetic 
Heisenberg model in the limit of large U: 

Hj = jJ2( s i- s j-m n j/*) ( 2 ) 

where J = At 2 /U. Si and rii are, respectively, the spin 
and the electron- number operators at the site i. Phys- 
ically speaking, the derivation of the Heisenberg model 
is as follows: The zeroth order effect of t for large U is 
to completely prevent the double occupancy of any site, 
thereby minimizing the large on-site interaction energy 
cost, which, in turn, gives rise to a low-energy Hilbcrt 
space in which sites are only singly occupied. This low- 
energy Hilbert space, however, is hugely degenerate since 
all states with single occupancy have exactly the same 
energy. Therefore, one has to investigate the next or- 
der effect of t. Exactly at half filling, there is no effect 
first order in t because, when it acts on the degener- 
ate Hilbert space mentioned above, the hopping term in 
Eq. © always creates a doubly occupied site which is out- 
side the low-energy Hilbert space. The antiferromagnetic 
Heisenberg model emerges through the second-order con- 
tribution which is the virtual hopping process (hence, 
J oc t 2 /U). Note that the virtual hopping process mini- 
mizes the kinetic energy cost in the presence of the strong 
on-site Coulomb repulsion. 

Away from half filling with addition of holes, the 
Heisenberg model generalizes to the t- J model: 

H t _j=V G (H t +Hj)V G , (3) 

where H j is given in Eq.© and 

ff t = -*5>lc 7 - (r + H.c.). (4) 

(i,3) 

The Gutzwiller projection operator, Vg, imposes the no- 
double-occupancy constraint, which handles the strong 
on-site interaction. It is important to note that, to or- 
der t 2 /U, the rigorous derivation of the t-J model from 
the Hubbard model requires the omission of three-site 
hopping termsi^, which can be justified at low doping 
because the hopping terms have zero matrix elements 
exactly at half filling and their contribution is roughly 
proportional to the hole concentration away from half 
filling. In some sense, there are two small expansion 
parameters involved in deriving the t-J model from the 
Hubbard model: t/U and the hole concentration, x. The 
three-site hopping terms may become sizable when either 
parameter is not small. 



As shown in the above, antiferromagnetism is a natu- 
ral consequence of the strong Coulomb repulsion at low 
doped regimes. This is, of course, consistent with exper- 
imental findings that there is a well-defined long-range 
antiferromagnetic order (also known as Neel order) at low 
doping. The next question, then, is when and if antifer- 
romagnetism can generate pairing; specifically, whether 
the t-J model contains superconductivity in realistic pa- 
rameter regimes (with non-zero doping). It is the goal 
of this paper to provide evidence for the existence for 
superconductivity in the t-J model. To this end, we in- 
vestigate the connection between the t-J model and the 
Gutzwiller-projected BCS Hamiltonian. Note that the 
Gutzwiller-projected BCS Hamiltonian is nothing but 
the BCS Hamiltonian in the presence of the strong on-site 
repulsion. 

The viewpoint advocated in this paper that antifer- 
romagnetism and high Tc superconductivity are inti- 
mately connected is shared by many of previous the- 
ories including the 50(5) theory 15 . The central hy- 
pothesis of the 50(5) theory is that antiferromagnetism 
and superconductivity are just two different manifesta- 
tions of a single object called superspin which combines 
the three-dimensional antiferromagnetic order parameter 
and the two-dimensional superconducting order parame- 
ter. Apart from the phcnomcnological appeal of this idea, 
it is crucial from the microscopic point of view that one 
should be able to derive an effective superspin model from 
well-known microscopic models such as the t—J model in 
a reliable manner. One of such attempts was to use some 
form of the renormalization-group transformation^. In 
this paper, however, we take a different approach, as out- 
lined below. 

This paper is organized as follows: The precise math- 
ematical form of the Gutzwiller-projected BCS Hamil- 
tonian is given in Sec. [H] where we discuss the physi- 
cal motivation for studying its connection to antiferro- 
magnetism. We also discuss the relationship between 
the ground state of the Gutzwiller-projected BCS Hamil- 
tonian and Anderson's resonating valence bond (RVB) 
stated (i.e., the Gutzwiller-projected BCS ground state). 
In Sec. IIIII by computing the wavefunction overlap via 
exact diagonalization in finite systems, we provide evi- 
dence suggesting that the ground state of the t-J model 
is closely connected to the ground state of the Gutzwiller- 
projected BCS Hamiltonian. In particular, we empha- 
size that, within the limits of numerical accuracy, the 
ground state of the antiferromagnetic Heisenberg model 
(i.e., the t-J model at half filling) is equivalent to that of 
the Gutzwiller-projected BCS Hamiltonian with strong 
pairing. In fact, this equivalence can be derived ana- 
lytically. In Sec. IIVI we provide an analytic derivation 
for the equivalence between the Heisenberg model and 
the Gutzwiller-projected BCS Hamiltonian with strong 
pairing. We conclude in Sec. 0by discussing physical 
implications of the Gutzwiller-projected BCS Hamilto- 
nian. 
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II. THE GUTZWILLER-PROJECTED BCS 
HAMILTONIAN 

We begin by providing a physical motivation for study- 
ing the connection between the t-J model and the pro- 
jected BCS Hamiltonian. For this purpose, it is conve- 
nient to write the pairing term of the BCS Hamiltonian 
in real space: 

-Hpair = C i1 C jl ~ C il C j'\ H.C. (5) 

which creates a resonance of the singlet pair between the 
site i and j: | tilj) — I W\j)- On the other hand, the sin- 
glet pair is energetically preferred by the antiferromag- 
netic exchange term in the t-J model: Si ■ Sj for J > 0. 
Therefore, despite their different appearances, the pair- 
ing term and the antifcrromagnetic exchange term seem 
to have a similar physical effect; they both prefer singlet 
pairs (at least, between nearest neighbors). This similar- 
ity serves as a motivation to investigate whether there is 
a connection between the t-J model and some form of the 
BCS Hamiltonian. In this paper, we take the Gutzwiller- 
projected BCS Hamiltonianii: 

-^bcs = 'PgHbcs'Pg 

= V G {H t +H A )V G , (6) 

where H t is given in Eq.Q. Also, the pairing term is 

H A = E A ^( c lT c ]l- c ac] T +H.c.) (7) 

(i.J-> 

where we are primarily interested in a pairing with d- 
wave symmetry, where Ajj = A if j = i + x and —A if 
j = i + y. However, we also examine the pairing with 
extended s-wave symmetry, where A !3 - = A for both 
j = i + x and i + y, and the Hamiltonian is denoted 
as H%cs- Note that Hq CS in Eq.@ has a very similar 
structure as H t .j in Eq.jSj). The only change is that Ha 
in Hq CS replaces Hj in H t -j, which is motivated by the 
possible correspondence between the pairing term and 
the antiferromagnetic exchange term. Note that H t + HA 
is the usual BCS Hamiltonian in the real space represen- 
tation, and hence Hbcs = H t + Ha- 

In order to provide evidence for the close connection 
between the t-J model and the Gutzwiller-projected BCS 
Hamiltonian, we will pursue two approaches with one 
being numerical (Sec. Illlf) and the other being analytical 
(Sec. IIVI) . But, before presenting our results, we would 
like to briefly discuss the relationship between the ground 
state of the Gutzwiller-projected BCS Hamiltonian and 
the RVB state proposed by Anderson^, which is just the 
Gutzwiller-projected BCS ground state: 

V'RVB = VnTg^bcs, (8) 

where Vg is the previously introduced Gutzwiller pro- 
jection operator imposing the no-double-occupancy con- 
straint. Vn is the electron-number projection operator 



which is necessary for constructing states with definite 
electron number from the BCS ground state, V'bcSj that 
has an intrinsic fluctuation in particle number. 

In general, V'rvb is not the exact ground state of H^ cs 
since the Gutzwiller projection does not commute with 
the BCS Hamiltonian, Hbcs: [V G ,H t + H A ] ^ 0. The 
difference between i/'RVB and the ground state of H^ cs 
is particularly severe at half filling because, as shown in 
the later sections, the ground state of H^ cs has long- 
range antiferromagnetic order (Neel order) at half filling, 
while ^rvb does not. Since the exact ground state of 
the t-J model has Neel order at half filling, it seems that 
V'RVB cannot be a good ansatz wavefunction for the t-J 
model, at least for regimes very close to half filling. This 
was one of the reasons why i/'rvb was proposed as a can- 
didate for the ground state of models with sufficiently 
strong "quantum frustration" (to destroy Neel order), 
examples of which include models with the next-nearest- 
neighbor exchange coupling, J', and antiferromagnetic 
models on the triangular lattice^. Also, f/'RVB was con- 
jectured to be a good ansatz wavefunction for the t-J 
model at moderate, non-zero doping. Its validity, how- 
ever, has remained very controversial even after many 
years of active research 18 i 19 i 20 i 21 i 22 i 23 i 24 . 

The situation is different for the ground state of H^ cs . 
In addition to the possession of Neel order at half filling, 
the ground state of Hq CS also has a very high overlap 
with the ground state of the t-J model at moderate dop- 
ing (in fact, unity overlap at half filling), as shown in the 
next section. Therefore, the ground state of H^ cs can be 
taken as a good ansatz wavefunction for the ground state 
of the t-J model at general doping. Now, this brings up 
an interesting question: how is the ground state of H^ cs 
related to the RVB state at non-zero doping? In Sec. 
it is argued that, despite severe differences at half fill- 
ing, the RVB state is, in fact, qualitatively similar to the 
ground state of H^ cs in doped regimes sufficiently away 
from half filling. 



III. NUMERICAL EVIDENCE 

In this section, we present numerical evidence for 
the close connection between the t-J model and the 
Gutzwiller-projected BCS Hamiltonian. While some nu- 
merical results are available in a previous paper by the 
author 17 , more details are given in this section along 
with new results. Before we provide detailed evidence, 
we would like to emphasize two important aspects of the 
numerical approach used in this paper. The first aspect 
pertains to the computing technique while the second is- 
sue concerns what to compute. 

First, we use exact diagonalization as the comput- 
ing technique. Exact diagonalization offers an unbi- 
ased approach, as opposed to other methods with bi- 
ased assumptions or guesses, which include the large 
N (or S) expansionMi 2 ^ and variational Monte Carlo 
simulationi&iiiSSiSii 2 ^. Uncontrolled approximations in 
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techniques mentioned above require an independent ver- 
ification of the involved assumptions. In contrast, ex- 
act diagonalization can determine the exact ground state 
wavefunction without bias 26 . The problem is, however, 
that its application is limited to studies of finite systems 
with relatively small spatial size. This limitation moti- 
vates the second aspect of the numerical approach in this 
paper. 

The second aspect of our numerical approach concerns 
what to compute. In general, evidence for long-range or- 
der is provided by relevant correlation functions. In the 
case of superconductivity, the relevant correlation func- 
tion is the pairing correlation function: 

F a0 (r - r') = (c\ (r)cj (r + a)q (r')c T (r' + 0)) , (9) 

where a, (3 — x,y. True off-diagonal long-range order 
(ODLRO) can be claimed only when F a p remains non- 
zero in the limit of large distance |r — r'|. Unfortunately, 
however, the small spatial size of finite systems accessible 
via exact diagonalization makes the distinction between 
true long-range order and short-range order (present even 
in normal states) ambiguous. Therefore, a measure of 
pairing order that is unambiguous even in finite system 
studies is needed. 

An inspiration comes from the fractional quantum Hall 
effect (FQHE), where our understanding of the subject is 
significantly advanced by the direct comparison between 
ansatz wavefunctions and exact states. It is well accepted 
by now that all essential aspects of the FQHE are ex- 
plained by the composite fermion (CF) theory 2 ! which 
is a general theory of the FQHE including the Laughlin 
stated as a subset. While there are various (both ex- 
perimental and theoretical) verifications of the compos- 
ite fermion theory, arguably the most significant is the 
amazing agreement between the exact ground state and 
the CF wavefunction: the overlap is practically unity for 
various short-range interactions including the Coulomb 
interaction 2 ^. Indeed, as was crucial in establishing the 
CF theory for the FQHE, we would in turn like to achieve 
the same methodological clarity for the t-J model. 

In order to give a perspective on the significance of 
wavefunction overlap in finite system studies, consider 
two randomly-chosen states in a Hilbert space with iVbasis 
basis states. Then, the possibility for having a large 
wavefunction overlap between those two states is roughly 
l/M^asis- So, roughly speaking, if the square of the over- 
lap between an ansatz state and the exact state is signif- 
icantly higher than 1/iVbasis) it can be argued that the 
ansatz state is a good representation of the exact state. 
The main finite system studied in this paper is the 4x4 
square lattice system, whose Hilbert space has 10 3 - 10 5 
basis states depending on the number of holes (even after 
translational symmetries are implemented as reported in 
this paper). Therefore, in our system, the possibility for 
having a large wavefunction overlap between two ran- 
dom states by chance is roughly 0.1% - 0.001%. We will 
show in the following sections that the square of the over- 
lap between the ground states of the t-J model and the 



Gutzwiller-projected BCS Hamiltonian is very close to 
unity at optimal parameter ranges (typically 90% - 100% 
depending on parameters). 

While it may seem straightforward at first to com- 
pute wavefunction overlap by exactly diagonalizing the 
t-J model and the Gutzwiller-projected BCS Hamilto- 
nian, there is a subtle, but physically important twist in 
applying exact diagonalization to the BCS Hamiltonian. 
The subtlety arises from the fact that there is a (coher- 
ent) fluctuation in particle number due to the pairing 
term in which has matrix elements mixing the Hilbert 
space of N e and N e ±2 electrons. In order to incorporate 
the particle-number fluctuation into finite system stud- 
ies, we diagonalize H^ cs in the combined Hilbert space 
of N c and N c — 2 electrons, which, in turn, invariably 
requires a careful treatment of the chemical potential. In 
essence, the chemical potential is adjusted so that the 
kinetic energy plus the chemical potential energy of the 
N e particle ground state is the same as that of the N e — 2 
particle ground state. Once the chemical potential is set 
this way, the mixing with other particle-number sectors 
such as the N e + 2 and N e ± 4 sectors can be shown to 
be negligibly small, even if it is allowed. For more de- 
tails, readers are referred to discussions in the following 
sections. 

Finally, the following notations are defined for future 
convenience: ipBCs(Nh, Nh + 2|7V) denotes the ground 
state of the Gutzwiller-projected BCS Hamiltonian ob- 
tained from the combined Hilbert space of Nh and Nh + 2 
holes in the system of N sites. (Note that the sum of 
the number of electrons and holes equals the number of 
sites: N e + N h — N). pN h =N denotes the number pro- 
jection operator which projects states onto the Hilbert 
space of states with Nq holes, and renormalizes the pro- 
jected states. ipt-j(Nh\N) is the exact ground state of the 
t-J model in the Hilbert space of Nh holes in N sites. 

Numerical evidence for the close connection between 
the ground states of the t-J model and the Gutzwiller- 
projected BCS Hamiltonian is presented as follows: In 
Sec. IIII Al the symmetry of the Gutzwiller-projected 
BCS Hamiltonian with d-wave pairing is discussed. 
Then, the wavefunction overlap between tpt-j{Nh\N) and 
the appropriately number-projected tp^ cs (Nh, Nh+2\N) 
is computed in Sec. IIII Bl for the case of undoped regime, 
which is followed by similar calculations in Sec. IIII CI for 
the optimally doped regime, and in Sec. IIII Dl for the 
over doped regime. 



A. Symmetry 

In general, before solving any Hamiltonian, one has to 
examine the symmetry of the Hamiltonian. Incorporat- 
ing symmetry is particularly important when two differ- 
ent Hamiltonians are compared since it is possible for 
their ground states to have completely different symme- 
tries from each other. The BCS Hamiltonian with d-wave 
pairing is particularly tricky in this respect because (i) 
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it does not conserve the particle number and (ii) it is 
not invariant under diagonal reflection, i.e., x <-» y (or, 
equivalently, under rotation in space by 7r/2). There- 
fore, it does not conserve parity with respect to diago- 
nal reflection. Surprisingly, however, the effects of the 
above two properties cancel and can be eliminated si- 
multaneously by applying the number projection opera- 
tor to the ground state of the Gutzwiller-projected BCS 
Hamiltonian, as explained below. (As far as parity is 
concerned, it is not important whether the BCS Hamil- 
tonian is Gutzwiller-projected; the Gutzwillcr projection 
commutes with diagonal reflection.) 

First, the non-conservation of particle number is, of 
course, a direct consequence of the pairing term in the 
BCS Hamiltonian, ejej + cj,cj. At first glance, the un- 
certainty in particle number may seem as a theoretical 
artifact in the BCS Hamiltonian. But, the coherent fluc- 
tuation in particle number is essential for superconduc- 
tivity since it is vital to the superfluid phase coherence. 
(The coherent fluctuation in superconductors stands in 
sharp contrast with the incoherent fluctuation in thermo- 
dynamic ensembles.) On the other hand, when a precise 
comparison with number eigenstates (such as the ground 
state of the t-J model) is required, the ground state of the 
BCS Hamiltonian should be projected onto the Hilbcrt 
space of states with a definite particle number. The num- 
ber projection operator, Vn, performs this task. 

Second, the non-conservation of parity with respect to 
diagonal reflection is due to the sign difference in <i-wave 
pairing amplitudes between the x and y direction. To 
be more precise, the non-conservation of parity is math- 
ematically expressed as follows: 

R d H£ cs (t, A)TZ d = Hg cs (t, -A) (10) 

where lZ d is the diagonal reflection operator and 
H^ cs (t, A) is the Gutzwiller-projected BCS Hamiltonian 
given in Eq.JSJ with the dependence on t and A shown 
explicitly Note the sign change of A in H^ cs in the right- 
hand side of Ea. ljTTl . Now, for convenience, Ea. lfTU|) is 
re-written as follows: 

K d Hg cs (t,A) = Hg cs (t, -A)K d (11) 

since 1Z\ — 1. What we want to prove in this section 
is that a number-projected eigenstate of H^ cs is also a 
parity eigenstate. 

We begin by defining \^f) as an eigenstate of 
#£cs(*>A): 

flf ca (t,A)|tf)=.E|*>. (12) 

Also, for convenience, let us write jvP) in vector form 
showing amplitudes in each individual particle-number 
sector. That is, 

(*l = (•■• ■■ ,c£ e+2 ,---), (13) 

where C N is the amplitude for the a-th basis element 
in the N e particle Hilbert space. Similarly, C N +2 is the 



amplitude for the /3-th basis element in the N e + 2 particle 
Hilbcrt space. 

Now, it is not too difficult to show that 

H$ ca (t,-A)\*)=E\V), (14) 

where E is the same energy as in Ea. (|12fl and 

<*! = (•■■, [-i]*ca., • •• . • • • ) ,(i5) 

where N e is restricted to even numbers since we are in- 
terested only in paired states. What accounts for the 
form of Eq. l|15l) is that the pairing term, Ha, always 
mixes the Hilbert spaces with particle numbers differ- 
ing by two. Attaching the relative sign according to the 
particle number (i.e., [— l]" 5 " and [— 1] 2 for the N e 
and N e + 2 Hilbert space, respectively) is tantamount to 
changing the sign of A in Ha- Note that a relative sign 
difference between distinct Hilbert spaces does not alter 
the hopping term since H t does not mix Hilbcrt spaces 
with different particle numbers. 

Now, let ltd act on both sides of Eq. i|12|) : 

it d Hg cs (t,A)\V} = Eil d \*), (16) 

which, with aids of Ea. i|ll|) . becomes 

Hi C3 {t,-A)K d \9) = ER d \iJf). (17) 

Then, by comparing Eq. IJ14J1 and l|17Jl , one is able to con- 
clude that 

ti d \V) = A|¥), (18) 

where A is a constant which is almost (though not quite 
yet) an eigenvalue. After acting on both sides by Vn=n c , 
the above equation becomes 

n d v N=Ne \^>) = xp n=n .\v). (19) 

(Note that T>n commutes with V, d - Also, there is no loss 
of generality in choosing N — N e .) Moreover, since 

Vn=nM = (-1)^^.1*), (20) 
Ea. (|19|) becomes 

K d V N=N J^) = X"P N=N M, (21) 

where A' = (— 1)~A. Consequently, Vn\^) is an eigen- 
state of the parity with respect to the diagonal reflection. 

Until now, only the spatial symmetry has been inves- 
tigated. (Note that H^ cs is invariant under the spa- 
tial translation, and therefore linear momentum is con- 
served.) Now, we would like to briefly consider spin ro- 
tation symmetry. The hopping term is obviously invari- 
ant under spin rotation. In addition, the pairing term 
can be also shown to be invariant under the spin rota- 
tion: [Ha, Stot] = where the total spin is S to t = J2i ^i- 
The essential physics of the spin-rotational invariance of 
the pairing term is the fact that Ha concerns only spin- 
singlet pairs, which are rotationally invariant. 
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B. Undoped regime (half filling) 

In this section, we provide numerical evidence in- 
dicating that, at half filling, the ground state of the 
Gutzwiller-projected BCS Hamiltonian is actually equiv- 
alent to the ground state of the t-J model in the limit of 
strong pairing, i.e., A/t — ► oo. To this end, the wavefunc- 
tion overlap between the two ground states is computed 
as a function of A/t via a modified Lanczos method for 
exact diagonalization. Note that the ground state of the 
t-J model is uniquely determined at half filling without 
any dependence on J/t, while the ground state of the 
Gutzwiller-projected BCS Hamiltonian depends on A/t. 
The reason for the former is that, at half filling, the t-J 
model becomes the antiferromagnetic Heisenberg model 
where J is just an overall scale factor. 

As mentioned previously, the main finite system stud- 
ied in this paper is the 4x4 square lattice system with pe- 
riodic boundary conditions. The 4x4 system is one of the 
most studied systems in numerical treatments^ because 
it is accessible via exact diagonalization, yet large enough 
to contain essential many-body correlations. We have 
checked that our results for the t-J model are in complete 
agreement with previous numerical studies^ for all avail- 
able cases. Using notations defined earlier in Sec. II I II the 
ground state of the t-J model is denoted as -0 t _j(O|16) 
for half filling: holes (16 electrons) in the 4x4 system. 
Also, the ground state of the Gutzwiller-projected BCS 
Hamiltonian at half filling is obtained by applying VN h =o 

to Vbcs(°> 2 I 16 )- Note that V'scst ' 2 I 16 ) is the g round 
state of the combined Hilbert space of and 2 holes (16 
and 14 electrons, respectively) in the 4x4 system. 

At half filling, the overlap between the ground states of 
the t-J model and the Gutzwiller-projected BCS Hamil- 
tonian is given by (Vw(0|16)| Av^ol^Icst ' 2| 16)> , the 
square of which is plotted in Fig. ^ as a function of A/t. 
As one can see from Fig. ^ t ne overlap is very high even 
for A/t ~ 1, and saturates very quickly to unity as A/t 
increases further. For sufficiently large values of A/t, 
the overlap is indistinguishable from unity to within nu- 
merical accuracy (We have actually studied A/t values 
ranging as high as 1000). Therefore, as far as our finite 
system is concerned, the ground state of the Gutzwiller- 
projected BCS Hamiltonian is exactly identical to the 
ground state of the t-J model in the limit of strong pair- 
ing, i.e., A/t — > oo. In view of the fact that there are 
roughly 10 3 basis states in the Hilbert space for the 4x4 
system at half filling even after translational symmetries 
are implemented, the high overlap is particularly salient. 
We have also checked that the above result also holds in 
the 4-site (2 x 2) and the 10-site (VlO x VlO) systems. It 
should be emphasized that the equivalence between the 
ground state of the t-J model and that of the projected 
BCS Hamiltonian with strong pairing does not necessar- 
ily mean strong superconductivity since, despite strong 
pairing, there is little charge fluctuation near half filling, 
and therefore little phase coherence. 

The equivalence between the ground states of the pro- 



(V-o(0|16)|^v„=oWg cs (0,2|16)) 




FIG. 1: Square of the overlap between the ground states of 
the t-J model at half filling (i.e., the 2D antiferromagnetic 
Heisenberg model) and the projected BCS Hamiltonian. Ex- 
act diagonalization is performed in the 4x4 system with 
holes (undoped regime). Inset: plotted is the square of the 
overlap in an extended range of A /t up to 10 showing quick 
saturation of the overlap to unity. 

jected BCS Hamiltonian and the t-J model at half filling 
has a very important physical implication for the origin 
of high Tq superconductivity supporting a long-standing 
conjecture^. The conjecture is that electrons are al- 
ready paired at half filling, and therefore form a conden- 
sate. However, electrons cannot superconduct (or, for 
that matter, even conduct) at half filling because there 
is no room for them to travel. But, away from half fill- 
ing, it seems natural that removing some fraction of elec- 
trons (i.e., doping) may trigger superconductivity by mo- 
bilizing electrons (and also causing charge fluctuations). 
In support of this idea, we show in Sec. 1111 CI that the 
overlap between the ground states of the projected BCS 
Hamiltonian and the t-J model remains very high even 
at non-zero doping for reasonable parameter ranges. 

C. Optimally doped regime 

In this section, we present numerical results for the 
wavefunction overlap between the ground states of the 
projected BCS Hamiltonian and the t-J model at a mod- 
erate, non-zero doping. Specifically, we study the situa- 
tion in which there are 2 holes in the 4x4 system, which 
roughly corresponds to the optimally doped regime. Us- 
ing notations defined earlier, the overlap is given by 

(■0t-j(2|16)|P JVh= 2|V'BCs( O ' 2 l 16 )) for the Hilbert space 
with 2 holes. An alternative representation of the 
overlap for this Hilbert space can be obtained from 
(ipt-j{2\16)\V Nh =2\ipBCs( 2 ^\ 16 )) which, however, leads 
to essentially the same conclusion, as shown later in this 
section. 

Fig. |21 displays square of the overlap defined by 
\(ipt-j(2\16)\V Nh = 2 \i>%cs( > 2|16))| 2 . Note that the over- 
lap now is a function of two parameters: A/t for the 
projected BCS Hamiltonian and J/t for the t-J model. 
For better visualization, a three-dimensional plot of the 
overlap is also given in Fig. [3J As shown in Fig. [21 and 
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FIG. 2: Square of the overlap between the ground states 
of the t-J model and the projected BCS Hamiltonian in the 
4x4 system with 2 holes (which roughly corresponds to the 
optimally doped regime). 




FIG. 3: Three-dimensional plot for the square of the overlap 
as a function of A/t and J/t in the 4x4 system with 2 holes. 
While this plot is essentially the same as Fig|5| the existence 
of two distinct regions of high overlap is more easily visualized. 



Fig. 13 there are two distinct regimes of high overlap: a 
weak-coupling regime (J/t < 0.08 and A/t < 0.1) and 
a strong-coupling regime (J/t > 0.08 and A/t > 0.1). 
These two regimes are qualitatively different in the sense 
that the symmetry of the ground states with respect to a 
spatial rotation by tt/2 changes from s-wave to d- wave at 
the boundary J/t ~ 0.08 for the ground state of the t-J 
model, and at the boundary A/t ~ 0.1 for the ground 
state of the projected BCS Hamiltonian. Note that, in a 
similar fashion for the case of parity with respect to di- 
agonal reflection, the angular momentum associated with 
spatial rotation by tt/2 can be shown to be conserved in 
the BCS Hamiltonian with ci-wave pairing. Note also 
that, because of the symmetry change, the overlap in the 
regime defined by J/t > 0.08 and A/t < 0.1 is precisely 
zero. 

The high overlap in the weak-coupling regime is rather 
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FIG. 4: Square of the overlap for the 4x4 system with 
2 holes, which is in the same regime as Fig|3 But, the 
ground state of the projected BCS Hamiltonian in this plot, 
"p JVh=2 |V'scs(2, 4|16)), is obtained from the combined Hilbert 
space of 2 and 4 holes in the 4x4 system, as opposed to the 
combined Hilbert space of and 2 holes in Fig|U It is impor- 
tant to note that qualitative features of this plot are basically 
identical to Fig. |2] while there are some quantitative differ- 
ences. 



trivial because, in this regime, both the projected BCS 
Hamiltonian and the t-J model Hamiltonian are reduced 
basically to the same Hamiltonian which is just the hop- 
ping Hamiltonian, H t , with the Gutzillwer projection; 
hence, we regard the equivalence at small J/t and A/t 
as a self-consistency check for the techniques used in this 
paper. The high overlap in the strong-coupling regime, 
on the other hand, cannot be trivially explained. As seen 
in Fig. [3 the maximum value of the overlap (as a func- 
tion of A/t) approaches unity as J/t increases, which can 
be also seen in Fig. [31 in the form of a rising ridge as J/t 
increases. The high overlap in the strong-coupling regime 
is therefore not accidental, but rather is connected to the 
unity overlap in the limit of strong coupling, which is 
in turn due to an intrinsic connection between antiferro- 
magnetism and superconductivity. 

Finally, as mentioned before, we study the overlap de- 
fined by (ip t .j(2\16)\P Nh=2 \^ cs (2, 4|16)). This alterna- 
tive definition is also valid because the ground state of the 
projected BCS Hamiltonian can be obtained in two differ- 
ent ways: one can apply Vn h =i either to the ground state 
of the combined Hilbert space of and 2 holes, or to that 
of the combined Hilbert space of 2 and 4 holes. While, 
strictly speaking, these two definitions become identical 
only in the thermodynamic limit, it would be an assur- 
ing self-consistency check of our approach if they produce 
similar results even in the finite system studies that we 
study. We find that this is indeed the case. Fig. ^depicts 
\(^ t .j(2\lQ)\V Nh=2 \i,% cs (2A\lQ))\ 2 showing that the es- 
sential features are basically identical to those in Fig. [21 
while there are some minor quantitative differences. 
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FIG. 5: Square of the overlap between the ground states of 
the t-J model and the projected BCS Hamiltonian in the 4x4 
system with 4 holes (which roughly corresponds to an over- 
doped regime). 



D. Overdoped regime 

We now present numerical results for the case of 4 holes 
in the 4x4 system, which roughly corresponds to the 
overdoped regime. In this regime, it is not strictly ap- 
propriate to regard the t-J model as the large-£7 limit 
of the Hubbard model due to the omission of three-site 
hopping terms mentioned in the introduction. Exper- 
imentally, however, superconductivity is weakened and 
eventually destroyed as the hole concentration increases. 
Hence, it is natural to expect that the overlap between 
the ground states of the projected BCS Hamiltonian and 
the t-J model becomes small in this regime. We show in 
Fig. [S] that this trend does indeed occur. 

Fig.0 displays K^j(4|16)]P A r h =4|Vg s(2.4|16))| 2 , 
which shows that the overlap is, in general, negligibly 
small except for the trivial case of J/t = and A/t = 0. 
Therefore, in the overdoped regime, the ground state of 
the projected BCS Hamiltonian is no longer a good rep- 
resentation of the ground state of the t-J model. 



IV. ANALYTIC DERIVATION OF THE 
EQUIVALENCE AT HALF FILLING 

It has been shown numerically in the previous section 
that, at half filling, the ground state of the Gutzwiller- 
projected BCS Hamiltonian is equivalent to the exact 
ground state of the t-J model in the limit of strong pair- 
ing, i.e., A/t — ► oo. While the numerical evidence is 
quite strong (the overlap between the two ground states 
is indistinguishable from unity to within numerical ac- 
curacy), questions regarding the validity of finite-system 
studies linger: especially, whether the overlap is actually 
equal to unity, or just very close to unity. If only the lat- 
ter is true, the overlap will diminish and will eventually 



vanish as the system size increases. 

We emphasize that, even if it is true that the over- 
lap does vanish in the thermodynamic limit, the high 
wavefunction overlap for the finite system still provides 
strong support for the ansatz wavefunction by demon- 
strating that it contains the correct physics. This sit- 
uation is, in fact, very similar to what happens in the 
FQHE: the overlap between the Laughlin wavefunction 
(the CF wavefunction, in general) and the exact ground 
state of the Coulomb interaction decreases as the system 
size increases, and eventually vanishes in the thermody- 
namic limit. The vanishing overlap in the thermody- 
namic limit is an inevitable consequence of the fact that 
the two ground states are not precisely identical. For 
example, while the Laughlin state at the lowest-Landau- 
level filling factor v — 1/3 is the exact ground state of the 
short-range interaction given by V 2 <5(r)2ii22, it is still an 
approximation for the exact ground state of the Coulomb 
interaction (relevant for experiments) albeit an extremely 
good one. 

In this paper, it will be shown that the overlap at 
half filling is actually unity in the strong-pairing limit; 
in other words, we will prove that, at half filling, the 
ground state of the Gutzwiller-projected BCS Hamilto- 
nian with strong pairing is identical to the ground state 
of the Heisenberg model. We will also prove a stronger 
statement that the two Hamiltonians do not merely share 
the same ground state, but also have in common the same 
low-energy physics. Note that the Hamiltonian for the 
Heisenberg model is Hj in Eq. J2J and the Gutzwiller- 
projected BCS Hamiltonian with strong pairing is given 
by i?B C g in Eq. © with t — 0, which is nothing but 



with the Gutzwiller projection: 

(i,3) 



,t J 



H.c.) V G . (22) 



Note that the strong-coupling limit, i.e. the large-A/i 
limit, is equivalent to the situation in which t — since 
tj A = is not a singular point. 



Gutzwiller projection as the large-f/ limit and 
Lieb's theorem 



We begin our analytic derivation for the equivalence at 
half filling by writing the BCS Hamiltonian with finite, 
repulsive on-site interaction U: 



H BCS +u(t, A ) = H t + H Ao + UJ2 



(23) 



where H t is defined in Eq. It is very important to 
note that H/\ is given in Eq. 10 with a bare pairing am- 
plitude, Ao, instead of a fully renormalized pairing am- 
plitude, A. The relationship between Ao and A should 
become clear at the end of the analytic derivation: we 
mention in advance, however, that A oc A\jU . 
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Our study of Hbcs+u is motivated by the fact that, 
in the large-t/ limit, Hbcs+u reduces to the Gutzwiller- 
projected BCS Hamiltonian. In addition, treating the 
Gutzwiller projection as the large-t/ limit instead of 
working directly with the projection facilitates the analy- 
sis. Motivated by numerical studies given in the previous 
section, we are particularly interested in Hrqs+u with 
strong pairing: 

Hao+u = Hbcs+u (t = 0, A ) 

- H Ao + Uj2 n n n U ( 24 ) 

i 

which, in the large-C/ limit, becomes in Eq. H22|) 
What we strive to prove is that, in the limit of large 
U/Aq, Ha +u bas exactly the same low-energy physics 
(i.e., the same ground state and the same low-energy exci- 
tations) as the Hubbard model, 7?Hub, m the limit of large 
U/t. This, in turn, means that the low-energy physics is 
identical for both H^ and Hj since iJnub becomes Hj 
in the laige-U/t limit. 

It is important to note a subtle, but very crucial differ- 
ence between the large-?/ behavior of the Hubbard model 
and that of the strong-pairing Gutzwiller-projected BCS 
Hamiltonian. To appreciate this distinction, let us re- 
write the Hamiltonian for the Hubbard model as follows: 

-ffHub = H t + U }^ riii n,ii 

i 

= H BCS +u(t,A = 0), (25) 

which shows that H^ub is a special case of Hbcs+u with 
A = 0, while Ha+u was also a special case of Hbcs+u 
but with t — 0. Hence, there is a parallel between iJjjub 
and Ha +u- Specifically, the equivalence between the 
low-energy physics of H^ and Hj can be recast as the 
equivalence between the low-energy physics of two differ- 
ent parameter points (i.e, between those of t/ Aq = and 
t/Ao = oo ) of the same Hamiltonian, Hbcs+u^t Aq), m 
the large-C/ limit. 

In fact, it has been shown by Affleck et al£& that 
two apparently very different mean-field solutions of the 
half- filled Hubbard model are actually equivalent. The 
first mean-field solution is considered by Baskaran, Zou, 
and Anderson^i as well as by Ruckenstein, Hirschfeld, 
and Appel2£ and later by Kotliai— , who performed a 
quadratic factorization of Hj assuming that (cifCji) ^ 
and (c\ a Cj a ) = for nearest neighbor i and j. On the 
other hand, the second mean-field solution, considered by 
Affleck and Marston^I, assumes the opposite situation: 
(°i1 c jl) = an d ( c la c jcr) 0- It thus seems that the 
mean-field solution with {c^Cji) ^= and (c^Cja) = 0, 
which corresponds to Hbcs+u with t = 0, is equivalent 
to that with (c^Cji) = and {c\ a Cj a ) ^ 0, which cor- 
responds to Hbcs+u with Aq = 0. In this paper, we 
prove that a stronger equivalence, not just between mean- 
field solutions, but also between the exact ground state 
of Hbcs+u with t = and with Aq = 0. We, however, 



emphasize that, despite their success in mean-field the- 
ory, approaches based on unitary transformations are not 
applicable in the exact treatment due to the fundamental 
difference in the large-C/ behavior of -ffnub and Ha +Ui 
which is discussed in greater detail below. 

Consider the effect of a large U in i?Hub- In this situ- 
ation, the zeroth-order effect of H t corresponds to keep- 
ing only the on-site repulsion term, in which case the 
ground state energy is exactly the same for arbitrary 
spin configurations as long as there is a single electron 
per site. Therefore, there is a huge 2 N degeneracy in the 
low-energy Hilbert space which is, in fact, the Gutzwiller- 
projected space. Now, let us investigate the next order ef- 
fect in the Gutzwiller-projected space. The first-order ef- 
fect of H t does not contribute to the Gutzwiller-projected 
space since, at half filling, the hopping term always cre- 
ates a doubly occupied site taking the state beyond the 
Gutzwiller-projected space. In other words, exactly at 
half filling, 

VcHtVc = 0. (26) 

The antiferromagnetic Heisenberg model, Hj, arises from 
the second-order virtual hopping processes. 

Next, we consider the effect of a large U in Ha +u- 
The zeroth-order effect of Ha is, of course, the same as 
that of H t in H-^nb'- the Gutzwiller-projected space is the 
low-energy Hilbert space. In contrast to -ffHub, however, 
there is a non-zero first-order effect of pairing term: 

Hi - V G H A V G , (27) 

where A is renormalized from the bare value Ao. Physi- 
cally speaking, the strong-pairing BCS Hamiltonian sets 
up the pairing resonance directly instead of by virtual 
processes. Therefore, it is not a priori clear why or if 
H^ has the same low-energy physics as Hj. In fact, the 
above real v.s. virtual contrast suggests that, if it exists, 
the equivalence between the low-energy physics of 
and Hj cannot be derived via a simple unitary trans- 
formation of the Hamiltonian. Instead, the equivalence 
must be connected to intricate physics. 

An essential point is that the large-C/ behavior of 
Hao+u a nd H^ub can be dealt with systematically in our 
approach where we analyze finitc-C/ models and take the 
large-C/ limit as the Gutzwiller projection. There is no 
singular behavior in the large-C/ limit of H Ao +u (because 
is non-zero). However, it is not obvious whether the 
ground state of the antiferromagnetic Heisenberg model 
is adiabatically connected to those of the Hubbard model 
with finite U. In fact, it can be proven by Nagaoka's 
theorem-2, that, infinitesimally away from half filling, the 
ground state of the Gutzwiller-projected hopping term, 
VcHtPd (which, in a naive sense, is the large-C/ limit 
of the Hubbard model) is ferromagnetic rather than an- 
tiferromagnetic. Therefore, one should be careful in how 
one takes this limit. Fortunately, there is a theorem by 
Lieb^i showing that the ground state of the Hubbard 
model is uniquely determined at any positive U and, as 
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a corollary, the ground state in the large-t/ limit is adi- 
abatically connected to those of finite U. Therefore, the 
ground state of the antiferromagnetic Heisenberg model 
is adiabatically connected to ground states of the Hub- 
bard model for which U is finite. (We assume that this 
connection is also valid for low-energy excitations.) 

We now present our analytic derivation step by step 
beginning with an outline summarizing each step: (i) We 
begin by separating Ha +u an d -^Hub into two parts: 
the saddle-point Hamiltonian (H.a +u and 7iHub, re- 
spectively) and the remaining Hamiltonian incorporating 
fluctuations around the saddle-point solution (SHa„+u 
and <57YHub, respectively). The saddle-point Hamiltonian 
is chosen so as to capture a possible singular effect of 
the spin-density-wave (SDW) instability at (tt,tt), which 
is inherent in Ha and H t in the presence of repulsive 
on-site interactions with arbitrary strength. 

(ii) The ground states of Ha q +u and 7^Hub are ob- 
tained for a general U. It is shown that, in both cases, 
the saddle-point ground state is separated from other ex- 
cited states by an energy gap proportional to U when 
U is large. It is also shown that these two saddle-point 
ground states become identical in the large-?/ limit. We 
denote this saddle-point ground state in the large-t/ limit 
as |Vo>- 

(iii) In the large-t/ limit, the low-energy Hilbert space 
for the full Hamiltonian (including the saddle-point and 
fluctuation part) is composed only of states that are con- 
nected to |-0o) via rigid spin rotations, {Hi}. We denote 
these states as {IV^)}- Then, we show that in the large- 
U limit, all matrix elements of SHa„+u and STinub are 
precisely the same in the low-energy Hilbert space with 
the same being true for the matrix elements of the two 
saddle-point Hamiltonians. 

(iv) Having shown that all matrix elements in the low- 
energy Hilbert space are precisely the same in the large- 
U limit, we argue that the strong-pairing Gutzwiller- 
projected BCS Hamiltonian and the 2D antiferromag- 
netic Heisenberg model have identical low-energy physics. 
The situation is depicted schematically in Fig. EJ 

Finally, the above discussion provides physical insight 
as to why the ground state of the Gutzwiller-projected 
BCS Hamiltonian is fundamentally different from Ander- 
son's RVB state at half filling. The difference originates 
from the singularity due to the SDW instability which oc- 
curs at half filling. This singularity cannot be captured 
by applying the Gutzwillcr projection to the BCS wave- 
function which, as constructed, intrinsically lacks long- 
range antiferromagnetic order. In fact, generally, no long- 
range order can be generated by applying local operators 
to states without long-range order. On the other hand, 
long-range antiferromagnetic order is fully captured in 
the ground state of the Gutzwiller-projected BCS Hamil- 
tonian, as shown in the analytic derivation. It has been 
shown in numerical studies that the semi-classical Neel 
configuration has a finite weight in the ground state of 
the Gutzwiller-projected BCS Hamiltonian 40 . Despite 
the difference at half filling, however, we argue in Sec. |V] 
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FIG. 6: Schematic diagram for the structure of the Hilbert 
space expanded by saddle-point basis states. The vertical 
separation between short lines denotes the energy difference 
between saddle-point states. A new sub-Hilbert-space S« is 
obtained by applying a rigid spin rotation Ri to every state 
in So- Without a loss of generality, one can choose So to be 
the Hilbert space with spin quantization axis aligned with the 
long-range spin order, tpi denotes the lowest-energy state of 
Si. Note that the above sub-Hilbert-spaces are linearly inde- 
pendent, but not orthogonal to each other (hence, overlapping 
ellipses). In the limit of large U, tpi's form an over-complete 
basis set for the Hilbert space at half filling, which is concep- 
tually very similar to the coherent-state basis set for a spin 
representation. 



that, sufficiently away from half filling, the ground state 
of the Gutzwiller-projected BCS Hamiltonian is, in fact, 
intimately related to the RVB state. 



B. Step (i) 

It is well known that the nesting property of the Fermi 
surface at half filling induces an instability toward long- 
range antiferromagnetic order (i.e., Neel order) in the 
ground state of the Hubbard model. Nesting is mathe- 
matically defined as follows: 

£k = -£k+Q, (28) 

where ek = — 2i(cos k x + cos k y ) is the kinetic energy due 
to hopping. Q = (it, tt) is the nesting vector. The strong- 
pairing BCS Hamiltonian with d-wave pairing symmetry 
has precisely the same nesting property for the gap func- 
tion: 

A k = — Ak+Q, (29) 

where Ak = 2Ao(cosfc :c — cosk y ) is the gap function for 
d-wave pairing. Therefore, in analogy with the Hubbard 
model, the nesting property of Ak generates a SDW in- 
stability at (tt, tt). Note that, for both Hamiltonians, the 
chemical potential is set to zero at half filling. Away from 
half filling, the chemical potential becomes non-zero, in 
which case the perfect nesting is ruined and the SDW 
instability disappears. 

In order to study spin orders, it is convenient to re- 
write the on-site interaction Hamiltonian in terms of spin 



11 



i i 

operator, S, = ^la^abCib- 



C. Step (ii) 



U^2nitn il = -^—^2sf + — ^2(n^+n i i) (30) 



It is convenient to re-write H-Aq+u m momentum 
space: 



where the last term becomes constant if the total number 
of electrons is fixed, and will henceforth be suppressed 
except when its consideration is necessary. 

Keeping in mind that there is an intrinsic SDW in- 
stability at half filling, we decompose the spin oper- 
ator into the stationary and fluctuation part: Si = 
(Si) + (Si — (Si)). Then, by retaining all terms up to first 
order in Sj — (Sj), one obtains the saddle-point Hamilto- 
nian, TLa +u- 

i 

- An fete], -etc!,. t +H.C.) 



(31) 



where the d-wave pairing symmetry is explicitly written 
in real space and (pi = — ^jf{Si)A +u is the spin expecta- 
tion value for the ground state of Ha +u- The remaining 
terms form the fluctuation Hamiltonian, 8TLa +u'- 



sn Ao+ u = -^E( s * + i^ 



(32) 



Similarly, the Hubbard Hamiltonian, i?Hub, can be de- 
composed into two parts: 

W H ub = -t^2(cl a c ja + R.c.) 

i i 
III / S \ 2 

SH Huh = -— y £ / {8 i +—<pA , (33) 



where <pi = — ^-(Sj)Hub 1S the s P m expectation value for 
the ground state of Wnub- 

Since we are interested in the SDW singularity at 
(tt,tt), we set both = 0n cos (Q • rj)z and ipi = 
ipo cos (Q • Yi)z with Q = (tt,it). Finally, the saddle- 
point solution is completed by the determination of the 
optimal value of 4>o (Po) by minimizing the ground state 
energy of Ha +u (^Hub) with respect to </>o (fo)- Note 
that the optimization of the ground energy amounts to 
the saddle-point condition for the auxiliary field in the 
path integral formulation^; hence the label "the saddle- 
point Hamiltonian" . 



T~Ia +u = 



d 2 k - 



(2tt) : 



^ c kT C -kl 



H.c. 



+ i/(^( c UA^h.,), 



(34) 



where the Einstein sum rule is applied on the a indices. 
Moreover, by defining a spinor field 'J'k, 



*w = 



( CkT \ 

C k+Q,T 
J 

-kl 



(35) 



V «U J 

one can re- write Ti-A a +u m a convenient 4x4 form: 



7~La +u 



8U 



d 2 k 



(27T) 



2 



where 



Mi 



^a x A k cr z 
A k a 2 &-a x 



(36) 



(37) 



The momentum integration, J n d 2 k, covers half of the 
Brillouin zone, which, for convenience, we choose to be 
the area in fc-space bounded by k y — k x ± 7r and k y — 
— k x ±n. 

While it is straightforward to obtain eigenvalues of 
Mk, it is instructive to diagonalize Mk in two steps. 
First, we apply the following unitary transformation onto 
-M k : 



Ml 



U 

o u- 1 



Mi 



u- 1 

U 



where 



U 



-( 1 1 



(38) 



(39) 



which also defines 7±(k): 



/ 7 +T (k) \ 
7_ T (k) 

-rli(-k) 



U 

o u- 1 



(40) 



Second, M'^ is diagonalized via Bogoliubov transforma- 
tion: 



a± T (k) 
4i(-k) 



«±(k) -«±(k)\ / 7±t(k) \, 
«±(k) u±(k) A7l;(-k)^ 4iJ 
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where 



where 



4(k)-4( k ) = 
-2«±(k)«±(k) = 



and 



£ k = JA^ + (0o/2)a. 



(42) 



(43) 



Then, it can be shown that the ground state is a vacuum 
of the Bogoliubov quasiparticles: 

a ± Am S A+u) = 0, (44) 
which is satisfied by the following wavefunction: 

\^a 0+ u) = n[ u +( k )+ v +( k )7i T ( k hii(- k )] 

ken 

x [ u _(k)+ u _(k) 7 l T (k) 7 i i (-k)]|0) .(45) 

Moreover, minimizing the ground state energy with re- 
spect to cj>Q leads to the saddle-point equation for </> : 



_3_ 

2U 



d 2 k 
(2n) 2 



1 



(0o/2) 2 



(46) 



The saddle-point equation in Eq. Ij46(l has a remarkable 
property that the solution exists for arbitrary values of 
U. Especially intriguing is that the solution is singular 
in the small-C/ limit: ln^o oc —Aq/U for U — > 0. As a 
consequence, there is long-range antiferromagnetic order 
(i.e., tfro 0) even for an arbitrarily weak interaction U , 
which is, in turn, adiabatically connected to the long- 
range order in the large- U limit. 

In the large-?/ limit, </>o is proportional to U: (j)Q = 
2C//3. The excitation spectrum E^, therefore, has a large 
energy gap proportional to U. Consequently, in the large- 
U limit, the ground state, I^Ao+c)' ^ s completely sepa- 
rated from other excitations as far as the saddle-point 
Hamiltonian, TLa +u, is concerned. Furthermore, one 
can show that the situation is exactly the same for the 
Hubbard model: the ground state of 7iHub is separated 
from other excitations by an energy gap proportional to 
U in the large-C/ limit. Low-energy fluctuations must, 
then, come from the fluctuation Hamiltonians, 6H.a +u 
and SHii u b- 

Before we discuss the effect of fluctuations, however, 
we analyze 7iHub to show that, in the large-C/ limit, the 
ground state of Ha q +u is identical to that of Hn u b- 
The saddle-point solution of the Hubbard model was 
discovered long ago 25 i 26 ; the ground state of Hnub is a 
fully filled Fermi sea state of the "negative-energy-mode 
quasiparticles" , associated with /?_ , while it is a vacuum 
for the "positive-energy-mode quasiparticles" , associated 
with f3 + . The precise definition of (3± is given by: 



//?+T(k) 
/3_ T (k) 



-m o 

m & oo 

f fc r] k 

V o o -% & 





f ^ \ 


) 


C k+Q,T 




«4 1 




\ Ck+Q,l / 



,(47) 



ik ~ Vk 



£k 



and 



£k = 



V / 4 + w^ 



The ground state of 7iHub is, then, given by 
= n^ T (k)/3l i (-k)|0). 



Kub) 



(48) 



(49) 



(50) 



kesi 



Similar to what was done in the case of the strong- 
pairing BCS Hamiltonian, minimizing the ground state 
energy leads to the saddle-point equation for tpo, which 
is given by 



3 

2U 



d 2 k 



1 



n (2^) 2 v/e k + (W2) 2 



(51) 



It can be shown from Eq. (|51|l that, in the large-?/ limit, 
ipo becomes proportional to U; more precisely (po = 2C//3, 
in which case £ fc = —r)k = 1/V% [See Eq. igSJ]. At the 
same time, for the strong-coupling BCS Hamiltonian, 4>o 
is also proportional to U in the large-C/ limit so that 
w+(k) = 1, u+(k) = 0, u_(k) = 0, and u_(k) = 1 [See 
Eq. (|42|l ]. which results in the following ground state: 



|^o) = n^ T (k)7l i (-k)|0). 



(52) 



ken 



As one can see, |V>o) is identical to the ground state of 
^Hub in Eq. I|50|) if the large-C/ limit is taken (Note that 
/3_(k) = 7-(k) in the large-C/ limit). 

It should not come as a surprise that (Si) is indepen- 
dent of parameters (such as U/t or U/Aq) since, in the 
large-C/ limit, the ground state, \ipo), itself is uniquely de- 
termined without the dependence on U/t or U/Aq. Re- 
member that, in the large-C/ limit, the Hubbard model 
becomes the Heisenberg model which has a single scale 
factor J. Therefore, the ground state (actually, any 
eigenstate) of the Heisenberg model is parameter free. 
The situation is similar for the strong-pairing Gutzillwcr- 
projected BCS Hamiltonian, H^. 

So, at least superficially, it seems that taking the large- 
U limit entails simply making U infinite and ignoring any 
effect of t or Ao- However, this is not correct because 
there is a reduction in the ground state energy due to 
finite A in Ha +u- The ground state energy of +u ) 
has the following form: 



N, 



= — K - 2 

8cr 



d 2 k 



A k + (^ /2)2 + ^,(53) 



where the last term C//6 comes from the constant term 
in Eq. (|30|l . Then, it can be shown that, in the large-C/ 
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limit. 



E 



Ag + U 



d 2 k ~ 



3A| 

U 
3A| 

U ' 



n ( 27 
d 2 fc 



A 2 

2 A k 



(2tt): 



■ (cos k x — cos /c a 



(54) 



where we use the fact that </>o = 2f7 /3 in the large-t/ limit. 
As one can see, for finite Ao, there is a reduction in the 
ground state energy, which is proportional to Aq/C/. As 
long as A is not zero, spin configurations are not random 
in the low-energy Hilbert space. This is the difference 
between the large-?/ limit and strictly infinite U . 

The situation is similar for the Hubbard model. The 
ground state energy of |"0H ub ) is given by: 



-^Hub 



7<Po 



2 2 



d 2 k 



N e 8U r " ~Ai(27r)' 
which, in the large-?/ limit, becomes 



(W2) 2 



U 



, (55) 



-^Hub 

N, 



d 2 k 2 

Jn WF^ 

3t 2 r d 2 k 



U Jn (2tt) 2 
3f 
' U ' 



(cos k x + cos k y Y 



(56) 



Therefore, again, spin configurations of the low-energy 
states arc not random. 



D. Step (iii) 

We have shown that the excitation spectra of Ha +u 
and Huub are gapped with a very large energy gap of 
0(U) in the large-?/ limit. However, the true low-energy 
excitation should be gapless at half filling, as required by 
Goldstone's theorem in two dimension where the spin 
rotation symmetry is broken by non-zero (S;). This 
dilemma is an artifact of the restriction of our attention 
to only a part of the full Hamiltonian. 

The saddle-point ground state, \ipo): is separated from 
other excitations of the saddle-point Hamiltonian. How- 
ever, when the full Hamiltonian (including the saddle- 
point and fluctuation part) is considered, there are many 
other states which have exactly the same energy as \ipo)'t 
these states are connected to |-0o) via rigid spin rotations. 
We denote these states as \tpi) = Ri|^>o) with {R^} being 
spin rotation operators. The energy of {l^i)} is the same 
for an arbitrary spin rotation because 

(il>i\H\1u) = (V'olRjtfRilV'o) = (MH\^o), (57) 

where H, either Ha+u or i?Hub, is invariant under spin 
rotation. In the large-?/ limit, therefore, the low-energy 



Hilbert space is composed only of Note that, 

although linearly independent, the basis states in {IV'i)} 
are not orthogonal to each other; forms an over- 

complete basis set, which is conceptually very similar to 
the coherent-state basis set for a spin representation^. 

Our task now is to evaluate matrix elements of the full 
Hamiltonian within the above low-energy Hilbert space. 
To this end, we first investigate the saddle-point equa- 
tion. Eq. (|51|l is very similar to Eq. (|46() . In fact, when 
all quantities are properly scaled, Eq. l|5T)l is exactly the 
same equation as Eq. i|46|) since the momentum, k, is 
a dummy integration parameter and Eq. I|51|l can be 
obtained from Eq. (|46|l by translating k by (it, 0). A 
constant translation in momentum does not affect the 
integral in the saddle-point equation. To be specific, if 
the solution of Eq. (|5T|l . <f>o/t, is a function of U/t: 



<p /t = f(U/t), 



(58) 



then the solution of Eq. I|46fl . 4>q/Aq, is related to U / Ao 
by the same function /: 

0o/A o = /(Z7/A ). (59) 

The precise functional form of f(x) is not important for 
general x, but it is useful to know that f(x) — 2x/3 when 
x > 1. 

Now, considering the definition of 5H.a +u and <57iHub, 
it is not too difficult to show that effects of both fluctua- 
tion Hamiltonians become equivalent in the large-?/ limit 
since (i) 4>q = (po for general U, if t is set equal to Ao, 
and (ii) low-energy Hilbert spaces for both models reduce 
to the same Hilbert space, in the large-?/ limit. 

So, in essence, it has been shown that, in the large-?/ 
limit, all matrix elements for both fluctuation Hamilto- 
nians are the same in the low-energy Hilbert space. It 
is very important to note that we are interested in the 
the lowest-order, non-zero effect of finite Ao and t in the 
limit of large U. Since, in the large-?/ limit, eigenstates 
themselves should be parameter free, lowest-order effects 
of finite finite Ao and t must emerge through their effects 
on and tpa, respectively. 

We now establish that matrix elements of both saddle- 
point Hamiltonians, Ha +u and Huub, are also identical. 
We begin by considering two arbitrary states from the 
low-energy Hilbert space of the BCS Hamiltonian with 
finite U, \rPj) = Rj\^ 0+U ) and |^.) = R k \^ 0+U )- The 
matrix element of Ti.A a +u between these two states is, 
then, 



3 



A 0+C /\^ n+ jR]Rck! r , 



A + U 



*}(E< 



A + U 
■ S,; ] Rfc 



^(vi; + a|{R]f^E^- s 4 



C 0+ A (60) 
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where Ef o+u = -3A 2 /U and \^ 0+u ) = |Vo) in the 
large-t/ limit. Note that the invariance of H\ Q and (j) 2 
under spin rotations is used. Similarly, in the large-?/ 
limit, the matrix element of 7iHub between the above 
states is as follows: 



^Hub(^|R]R^0 

ipo 



o\ ^0 



(61) 



where E^ uh = —3t 2 /U in the large-C/ limit. Therefore, it 
is clear from a comparison of Eq. (|60|l and l|61l) that the 
matrix elements are identical, provided that the large- 
U limit is taken while t is set equal to An. Note that 
<f>i — ifi for general U if t = An. 



E. Step (iv) 

We have learned two lessons from Step (ii) and (iii). 
First, the ground states of the saddle-point Hamiltoni- 
ans, HHub and H.a +u, are identical in the large-?/ limit, 
and therefore low-energy Hilbert spaces (which are ob- 
tained by applying rigid spin rotations to the saddle- 
point ground states) are also identical in this limit. Sec- 
ond, in the low-energy Hilbert space, matrix elements 
for the fluctuation Hamiltonians are identical for both 
the Hubbard model and the strong-pairing BCS Hamil- 
tonian with the same being true for the two saddle-point 
Hamiltonians. Therefore, we conclude that the strong- 
pairing Gutzwiller-projected BCS Hamiltonian at half 
filling is equivalent to the 2D antifcrromagnctic Heisen- 
berg model. 



Relation to the projected BCS Hamiltonian 
with extended s-wave pairing: a corollary 
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FIG. 7: Square of the overlap between the ground state 
of the antiferromagnetic Heisenberg model, and those of 
the Gutzwiller-projected BCS Hamiltonian with Ak = 
2 A (cos k x +cos k y ) [extended s-wave pairing] and 2 A (cos k x — 
cos [d-wave pairing]. 



Gutzwiller-projected BCS Hamiltonian with extended s- 
wave pairing is also equivalent to the antiferromagnetic 
Heisenberg model. 

To test this assertion numerically, we compute the 
wavefunction overlap between the ground states of the 
Heisenberg model and the strong-pairing Gutzwiller- 
projected BCS Hamiltonian with extended s-wave pair- 
ing. As before, we analyze the 4x4 square lattice sys- 
tem via exact diagonalization. To be specific, the strong- 
pairing Gutzwiller-projected BCS Hamiltonian with ex- 
tended s-wave pairing is given by: 



n sBCS 



V G (H t + H sA )V G , 



where 



A Vfctct . . - etc!. . . 

/ A i\ i+x,l %l i+x,T 



+3/4 »1 *+y,T 



(62) 



H.c.) 

H.c.). (63) 



In this section, we would like to discuss an interest- 
ing corollary of the analytic derivation given in previ- 
ous sections. Our derivation is actually valid also for 
the strong-pairing BCS Hamiltonian with an extended 
s-wave pairing: A k /A = 2(cosfc x + cosk v ). This is due 
to the fact that the saddle-point equation for extended 
s-wave pairing is the same as that of d-wave pairing and, 
consequently, all matrix elements remain the same as ob- 
tained for d-wave case. Note that extended s-wave pair- 
ing is induced by nearest-neighbor pairing amplitudes 
with the same sign for the x and y directions in contrast 
to the opposite signs in d-wave pairing. However, despite 
this important distinction, it is proven as a corollary to 
our d-wave result that, at half filling, the strong-pairing 



Fig. shows that, for sufficiently large A/t, the over- 
lap indeed becomes essentially unity for both Ak = 
2A(cosfc a; + cosky) [extended s-wave pairing] and 
2A(cos k x — cos ky) [d- wave pairing]. 

As before in the d-wave case, it is also interesting to 
examine the doped regimes. While there is no difference 
between extended s-wave and d-wave pairing at half fill- 
ing, the t-J model prefers d-wave pairing at finite doping 
(at least in the k = (0,0) channel). Fig. [S] shows the 
wavefunction overlap for the 4x4 system with 2 holes. 
Except for the trivial parameter regime of very small J/t 
(top panel in Fig. [SJ, the overlap essentially vanishes for 
extended s-wave pairing in contrast to the high overlap 
for the case of d-wave pairing (bottom panel in Fig. |8J. 
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FIG. 8: Square of the overlap for the 2-hole case. Plotted is 
the square of the overlap between the ground state of the t — J 
model and those of the Gutzwiller-projected BCS Hamilto- 
nian with Ak = 2 A (cos k x + cosfc H ) [extended s-wave pairing] 
and 2A(cos k x — cos k y ) [d-wave pairing]. 



V. CONCLUSION 

We conclude by discussing the physical implications 
of the close connection between the Gutzwiller-projected 
BCS Hamiltonian and the t-J model. There are prac- 
tical implications as well as fundamental ones. Funda- 
mentally, the close connection between the Gutzwiller- 
projected BCS Hamiltonian and the t-J model provides 
evidence for the existence of superconductivity in the t-J 
model, which, in turn, may suggest that it is not only 
possible, but rather natural that the pairing in high Tc 
superconductors is caused purely by electronic correla- 
tions. 

Now that we have reason to believe that the anti- 
ferromagnetic interaction is closely related to electron- 
electron pairing, the next natural question is how 
this connection can be used for the quantitative un- 
derstanding of experiments. While it is in principle 
straightforward to apply exact diagonalization to the 
Gutzwiller-projected BCS Hamiltonian, exactly solving 
the Gutzwiller-projected BCS Hamiltonian is as compli- 
cated as solving the t-J model in the first place. The 
true practical advantage of analyzing the Gutzwiller- 
projected BCS Hamiltonian, however, lies in the fact that 
the Gutzwiller-projected BCS Hamiltonian is much more 
amenable to the ansatz wavefunction approach than the 
t-J model. 

Ansatz wavefunction approaches have been very suc- 
cessful in various strongly correlated problems in con- 
densed matter physics. Two very salient examples are 
the fractional quantum Hall effect (FQHE) and the liquid 
Helium. Despite their apparent difference in the physical 
context, there is actually a rather profound commonality 
between the above, two examples. In both cases, dom- 
inant short-range correlations are first separated from 
long-range correlations. Then, effects of the strong short- 
range correlations are captured by the attachment of the 



Jastrow factor. The specific functional form of the Jas- 
trow factor, of course, depends on the nature of the prob- 
lem at hand. After short-range correlations are taken into 
account via the Jastrow factor, it is assumed that resid- 
ual long-range correlations are relatively much weaker so 
that the rest of the ansatz wavefunction describes essen- 
tially non-interacting quasiparticles. 

Specifically, the ansatz wavefunction for the FQHE is 
given by the composite fermion (CF) wavefunction 27 . 
In essence, the CF wavefunction is a product of the 
Jastrow factor, ni<j( z « — z j) 2p with p an integer and 
Zj = Xj + iyj , and a non-interacting fermionic wavefunc- 
tion of new quasiparticles known as composite fermions. 
It is known that the above functional form of the Jas- 
trow factor is most effective in minimizing the Coulomb 
repulsion between electrons in the lowest Landau level; 
the wavefunction overlap between the CF wavefunction 
and the exact ground state is practically unity for all 
available finite-system studies^. 

Another example is liquid Helium. The microscopic 
wavefunction for normal state 3 He is written in terms of 
an algebraic product of the Jastrow factor and the Slater 
determinant of plane- wave states 41 -: 



exp 



i<j 



u{ri 



Dot 



(64) 



where u(r) is determined by the interaction between He 
atoms: typically, u(r) cx 1/r 5 . While the Slater determi- 
nant provides the necessary antisymmetrization as well 
as the low-energy physics, the Jastrow factor captures 
dominant short-range correlations. Similarly, the ansatz 
wavefunction for superfluid 4 He is given by: 



4 Ho 



exp 



(65) 



l<3 



where u'(r) is basically identical u(r) except for slight 
changes due to the mass difference between 3 He and 4 Hc. 
Since the ground state wavefunction for non-interacting 
bosons is just a constant, the ansatz wavefunction for su- 
perfluid 4 He can be also viewed as the product of the 
Jastrow factor and a wavefunction for non-interacting 
bosons. 

From the discussion so far, it seems promising that a 
range of strongly correlated problems can be attacked 
via the Jastrow-factor approach, provided that resid- 
ual, long-range correlations are weak so that, to a good 
approximation, they can be treated as those of non- 
interacting quasi-particles. One criterion for the weak- 
ness of residual interaction may be whether or not the 
residual interaction can cause an instability toward a new 
phase. In other words, as long as the ground state for the 
original interaction (causing the strong short-range cor- 
relation) is not completely different from that of weakly- 
interacting (emergent) quasi-particles, the Jastrow-factor 
approach can provide a good ansatz wavefunction. 

The situation for the t-J model, or the Hubbard model 
in the large- J7 limit, is quite intriguing. First, note that 
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FIG. 9: Square of the overlap for the case of 2 holes in the 
y/W x vTO system. Plotted is the square of the overlap 
between the ground state of the t — J model and those of 
the Gutzwiller-projected BCS Hamiltonian with 2A(cos k x — 
cos k y ) [d-wave pairing]. The discontinuity in the square of the 
overlap is due to the sudden symmetry change in the ground 
state. 



the Gutzwiller projection can be actually regarded as the 
Jastrow factor imposing the no-double-occupancy con- 
straint that occurs because of the strong on-site repul- 
sion in the large-t7 limit. Keeping in mind that the t-J 
model is the Hubbard model in the large-C/ limit, it may 
be conjectured that the ground state of the t-J model 
is described well by the Gutzwiller-projected Fermi-sea 
state. Unfortunately, however, this is not true at half 
filling. We have learned in Sec. IIV CI that, at half filling, 
even the infinitcsimally weak interaction, U, can cause 
a spin-density-wave (SDW) instability leading to long- 
range antiferromagnetic order, i.e., Neel order. There- 
fore, at least close to half filling where Neel order per- 
sists, the Gutzwiller-projected Fermi-sea state cannot be 
a good ansatz wavefunction. 

On the other hand, away from half filling, the SDW in- 
stability disappears because perfect nesting is ruined with 
non-zero doping, as shown in Sec. IIVBI One may assume 
now that the Gutzwiller-projected Fermi-sea state can be 
a good ansatz wavefunction for finite, non-zero doping. 
This is, however, not true either because there is another 
instability caused by electron-electron pairing. It is, in 
some sense, the goal of this paper to show that the t-J 
model contains a pairing instability. Motivated by the 
fact that the Gutzwiller projection can play a role of the 
Jastrow factor for the Gutzwiller-projected BCS Hamil- 
tonian, we now show that the Gutzwiller-projected BCS 
wavefunction (i.e., the RVB state) is, in fact, a good 
ansatz wavefunction for the Gutzwiller-projected BCS 
Hamiltonian at finite, non-zero doping, which, in turn, 
leads to the eventual connection between the ground 
state of the t — J model and the projected BCS wave- 
function. 

While there have been studies showing that, for fi- 
nite doping, the Gutzwiller-projected BCS wavefunction 
yields good agreement with experiments^ as well as with 
some numerical studies^, we provide a more direct piece 
of evidence. Hasegawa and Poilblanc-S have shown that 



half filling 


"4>tJ = ^BCS / ^RVB 


moderate doping 


t/jtJ ~ ^BCS — ^RVB 



TABLE I: Interrelationship between the ground state of the 
t — J model, tptj, the ground state of the Gutzwiller-projected 
BCS Hamiltonian, V'bcsi an d the RVB state, V'rvb- 



the overlap between the projected BCS wavefunction and 
the ground state of the t — J model is high (~ 90%) for 
the case of 2 holes in the -s/lO x VTO lattice system. In or- 
der to make a direct comparison, we have computed the 
overlap between the ground state of the projected BCS 
Hamiltonian and that of the t — J model in the same 
system; we find that the optimal overlap for the ground 
state of the projected BCS Hamiltonian is more than 
98% for the same y/lO x VTO lattice system, as shown in 
Fig. El Therefore, it is shown (at least, in finite system 
studies) that the projected BCS wavefunction is a good 
ansatz wavefunction for the projected BCS Hamiltonian, 
as expected from the Jastrow-factor approach. 

Regarding the work of Hasegawa and Poilblanc 19 , it is 
interesting to note that they also computed the overlap 
between the projected BCS wavefunction and the ground 
state of the t — J model in the 4x4 lattice system and 
they found that it was very low (~ 4%). This sudden 
drop in the overlap is a finite-system-size artifact due to 
the fact that the half-filled Fermi sea state (technically 
speaking, its chemical potential) is ill-defined in the LxL 
system with L an even integer (e.g., the 4x4 system) , 
while there is no such problem in the L x L + 1 system 
with L an odd integer (e.g., the VTO x y/W system when 
L — 3). It is important to note that the ground state of 
the projected BCS Hamiltonian in our approach does not 
suffer from the above problem regardless of site numbers. 

The close relationship between the projected BCS 
wavefunction and the ground state of the projected BCS 
Hamiltonian at moderate doping is rather important 
since it provides additional support for the existence 
of superconductivity in the projected BCS Hamiltonian 
and, eventually, for that of the superconductivity in the 
t — J model through the connection between the pro- 
jected BCS Hamiltonian and the t — J model; large-scale 
Monte Carlo simulation studies have shownSLSI that the 
projected BCS wavefunction has the long-range pairing 
correlation in the thermodynamic limit. In this perspec- 
tive, the goal of our work in this paper is to show why 
the projected BCS wavefunction can be a good ansatz 
wavefunction for the t — J model at moderate doping, 
while it is not true at half filling. The interrelationship 
between the ground state of the t — J model, that of 
the Gutzwiller-projected BCS Hamiltonian, and the RVB 
state is schematically shown in Table Q] 

It is interesting to note that the recently proposed state 
for the Gossamer superconductor* 4 ^ 3 , is actually the BCS 
wavefunction with a partial Gutzwiller projection allow- 
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ing the double occupancy with suppression weights: 
^partial = JJ(1 - an iT n f .|.) 

i 

= e-^f^fHi, (66) 

where = 1 — a. The above expression suggests a 
very natural interpretation that the partial Gutzwiller 
projection is, in fact, also a specific form of the usual 
Jastrow factor. Incidentally, the Fermi sea state with 
the partial Gutzwiller projection, -pP artlal ; was studied 
previously in the context of the metal- insulator transition 
in the two-dimensional Hubbard models*. 

Finally, it is also interesting to note that, in one di- 
mension where there is no instability for either spin- 
density-wave or pairing, the exact ground state of 
the antifcrromagnctic Heisenberg model (i.e., Bethe 



ansatz solution 45 ) is actually very closely related to the 
Gutzwiller-projected Fermi-sea stateiS*^^, as expected 
from the Jastrow- factor argument. 
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